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Abstract 

This paper presents a guidance concept 
for use on the first stage of a candidate Mars 
Ascent Vehicle (MAV). The guidance is based 
on a calculus of variations approach similar to 
that used for the final phase of the Apollo 
Earth return guidance. A three degree-of- 
freedom (3DOF) Monte Carlo simulation is 
used to evaluate performance and robustness 
of the algorithm. 

1 Nomenclature 

£ Energy ( m 2 /s 2 ) 

g Gravitational acceleration (m./.s 2 ) 

H Hamiltonian 

% Angular Momentum (m 2 / s) 

h Altitude (m) 

L Lift (AO 

m Vehicle mass (kg) 

R e Planetary radius (m) 

r Radius from planet center (m) 

T Thrust (N) 

t Time (s) 

V Velocity (m/s) 

a Angle of attack (Thrust Angle) (deg) 

(3 Out-of-plane thrust angle (deg) 

7 Flight path angle (deg) 

S Variation 

d Partial Derivative 

A Costate 

// Gravitational parameter ( m 3 /s 2 ) 

v Latitude ( deg) 


t Dummy variable for time 

<j) Cost functional 

'(/’ Heading angle (deg) 

BO Value at Burnout 

a Value at Apoapsis 

targ Target Value 


2 Introduction 

Although efforts are underway to con- 
tinue improvements in reliability and sensi- 
tivity of robotic planetary probes, they will 
not, in the forseeable future, be able to 
match the examination and analysis capabil- 
ities available here on Earth. One solution 
to this dilemna is to retrieve planetary sam- 
ples for analysis here. This has been proposed 
for samples from Mars starting with the 2003 
launch opportunity. 

One plan calls for a lander to be sent 
to Mars to collect soil samples and launch 
them into orbit around Mars. The samples 
will remain in orbit until the orbiter and lan- 
der launched in 2005 reach Mars. The 2005 
lander will collect more samples and launch 
them into Mars orbit, and the 2005 orbiter 
will then rendezvous with both sets of sam- 
ples and return them to Earth. 

The Mars Ascent Vehicle (MAV) is sub- 
ject to severe design constraints. In ad- 
dition to the usual premiums on weight, 
volume, and budget, the MAV must oper- 
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Figure 1 . MAY Configuration. 


ate somewhat autonomously after being sub- 
jected, unattended, to a severe environment 
for nearly a year. As a result of these and 
other constraints, The MAV has unique chal- 
lenges in its design, especially for guidance 
and control. 

Because of the long travel times, the 
MAV will have solid-fuel motors. Figure 1 
shows the general configuration considered in 
this report. 

The ascent trajectory from the Martian 
surface begins with a high thrust phase that 
lasts approximately 20 seconds. The MAV 
then coasts for approximately 200 seconds, 
at which time it repoints, spins up to 20 rpm, 
separates the spent first stage and fires the 
second stage. One half orbit later, the third 
stage motor (which is mounted backwards to 
the other stages) is fired to circularize the fi- 
nal orbit. The second and third stages are not 
guided, though the repointing maneuver may 
be modified to account for an off-nominal first 
stage burn/coast. This sequence is illustrated 
in figure 2. 

While the first stage motor is burning, 
the vehicle is controlled by vanes in the rocket 
exhaust. After first stage burnout, aerody- 
namic surfaces are available to reorient the 
vehicle, but because of the low density of 
the Martian atmosphere, the ability to ad- 
just the first stage trajectory is limited. Thus, 
the objective of the first stage guidance is to 
achieve the highest degree of accuracy in the 
desired burnout conditions, subject to uncer- 
tainties in the winds, atmospheric density, ve- 
hicle/payload mass, and total impulse of the 
motor. The short burn time requires that the 


Reorient, 
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stage 2 burn 



1/2 orbit coast 


Figure 2. MAV Ascent Sequence. 

first stage guidance be very fast and robust 
to a rapidly changing plant. 

The scheme employed for guidance dur- 
ing the first stage uses an approach similar 
to that used for the final phase of the Apollo 
Entry Guidance [1, 2, 3]. Based on a nom- 
inal trajectory, the sensitivities of the final 
state (here the burnout state) to variations 
in the current state are determined and used 
to drive those variations to zero at the final 
time. 

The next section gives the derivation of 
the feedback equations for both the in-plane 
and out-of-plane components. Section 4 de- 
scribes the implementation of the algorithm 
in a numerical simulation and describes some 
results of that implementation. 

3 Theoretical 
Development 

The in-plane differential 

dynamical equations for a rocket ascent are 
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as follows. For altitude, 

h = V sin (7) 


( 1 ) 


• Tcos(a) . 

V = 0Sin(7) 


m 


and the equation for flight path angle is 

V 2 x 


= rsin(q) 
; mV 


R e + h 


H 


\ h V sin(7) + Ay 


T cos(cr) 


—Ay g sin (7) + A 7 
V 2 cos (7) 

7 (Re + h) 

and the costate equations are: 

A h = 


m 

T sin (cc) 
mV 

- X^g cos (7) 


dH_ 

dh 


A 7 D 2 cos (7) 
(Re + h ) 2 


A 


V 


dH 

~dV 

—Aft sin (7) - 

A 7 T sin (a) 
+ mV 2 


2A 7 D cos(7) 
R e T h 


\ — 


dH 

c?7 

-\ h V cos (7) + g\ v cos (7) 
A 7 V" 2 sin(7) 


R e + h 


- X 1 g sin (7) 


From [1], assuming that the perturbation in 
the control will be constant, 


where h is the altitude, V is the velocity, and 
7 is the flight path angle. The equation for 
velocity is 


Sa 


—A T (t)6x(t) 


A a (t) 


where 


(2) 


A, 


«(*) = - f f xT ( t ) 


df(r) 


dr 


or 


- g cos (7) (3) 


A a (t) = 


da(r) 

A yT sin (a) A 7 T cos (a) 


m 


mV 


(8) 


(9) 


( 10 ) 


The function to be maximized is the energy 
at burnout, so the Hamiltonian is [4]: 


(4) 


The above equations (5), (6), (7), (10) 
are integrated backwards from the final con- 
dition using states from the nominal trajec- 
tory. It is desired to match apoapsis to put 
second stage burn at the right position. So, 
let 

0 (X a I'targ) (H) 

The boundary conditions for the costates are: 


(5) 


A y(tB0) ~ 
A h(tBo) = 


d<j> 

dj 

d<p 

dh 


— 2(r a Ttarg) 


t=t B o 


2 (fa r targ ) 


t=tBO 


dr a 

dj 

dr a 
dh 


t—tso 

(12) 


A (+ 1 _ d( t> 


2 (r a - rtarg 


t t BO 

(13) 


(6) 


t=t B o 

From orbital mechanics: 

Veo^bo cos(7bo) = V a r a 


and 


V 2 

v BO 


h 


V a h 


t=t B o 

(14) 


(15) 


(16) 


(7) 


2 r B o 2 r a 

where the subscript 0 denotes quantities at 
burnout and the subscript a denotes quanti- 
ties at apoapsis. Define the energy and angu- 
lar momentum 
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V 2 

v BO 


h 

r B o 


(17) 


3 



'H — V BO^ BO COs(t bo) (18) 


Solving, 


m 

dj 


—r bcA BO sin Ibo 


(30) 


T BoVbO COS 7 bo 


n 

K 


(19) 


d£ 



( 31 ) 


substitute eqn(19) into eqn(16) and solve 
for V a : 


V a 


n±Vi £ 2 + 2 sn 2 

n 


( 20 ) 


The higher velocity is the periapsis root. 
The minus root is desired. Note: Energy will 
be negative for an elliptic orbit. The velocity 
at apoapsis will still be positive. 


V a = 


H~ vV + 2£U 2 

n 


( 21 ) 


Substituting into eqn(19) and simplifying 

n 2 


r n = 


H - vV 2 + 2£U 2 


(22) 


Note that the final state will never de- 
pend on the final control so the “control 
costate”, A a , will always have a final condi- 
tion of 0. When this is used in equation (8) it 
implies that a state perturbation at the final 
time requires infinite control to be corrected 
at the final time, i.e. instantaneously. 

Also note that, for flight implementa- 
tion, very little of this process occurs on- 
board. The costates, A 7 , Ay, etc, are de- 
termined from a nominal trajectory prior to 
flight. The costates are stored as tables or 
polynomials and the control is then a simple 
function of these stored values and the cur- 
rent state. 

The out-of-plane equation is 


The boundary conditions on the costates are 
then: 


dr a 1 d'H U dV a 

dh ~ V a dh V a 2 dh 

(23) 

on 

7 ^ — VBO cos 7 BO 

(24) 

d£ n 

dh r 2 BO 

The Vbo and 7 bo derivatives are 
analogous, with 

(25) 

exactly 

dr a 1 d'H n dV a 

dV ~ V a dV V 2 dV 

(26) 

d'H 

7^7 = f BO COS 7 BO 

(27) 

d£ -V 

aV - Vbo 

(28) 

dr a 1 d'H H dV a 

d^i v a dj V 2 97 

(29) 


• _ Lsm(<p) + T sin(/3) 
mV cos (7) 

V 

cos (7) cos(0) tan(z/) 

where v is the latitude, (3 is the out-of-plane 
thrust angle (similar to sideslip angle), and 
L is the lift force. It is assumed that the lift 
force is negligible compared to the T sin (f3) 
term and that the entire flight takes place 
near enough to the equator that the last term 
can be neglected. Writing as a finite differ- 
ence and solving for sin (/?), results in 

smffl = mVc ° S ^ (33) 

where ip n ew is the commanded heading angle 
and 'll) M is the current heading angle. For this 
implementation, the command was chosen as 
a ramp in time that brings the nominal tra- 
jectory to a 45 degree inclination at burnout. 
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4 Numerical Simulation 
Results 

The guidance algorithm described above 
was implemented in a numerical 3DOF simu- 
lation using the Program to Optimize Simu- 
lated Trajectories (POST) program [5]. The 
simulation was incorporated into a Monte 
Carlo analysis with dispersions as listed in 
Table 1. 


Table 1. Monte Carlo Inputs 

Variable 

Range 

Dist 

Launch Altitude 

0-2 km 

U 

Launch Latitude 

±0.1 deg 

G 

Launch Longitude 

±0.1 deg 

G 

Launch Azimuth 

41.74 ±1.8 deg 

G 

Launch FPA 

48.91 ±1.8 deg 

G 

E-W Wind 

±50 m/s 

U 

N-S Wind 

5-30 m/s 

U 

Propellant Mass 

38.418 kg ±0.3% 

G 

Payload Mass 

2.80 ±0.4 kg 

G 

Thrust 

5872.0 N ±4.0% 

G 

lsp 

279s ±1% 

G 

C A 

±5% 

G 

Cn 

±5% 

G 


The first column of Table 1 lists the 
quantities that were dispersed within the lim- 
its shown in the second column. The final 
column denotes the type of random distri- 
bution sampled; ’G’ for Gaussian and ’U’ 
for uniform. Random atmosphere variations 
were also included based on MarsGRAM [6]. 
The simulation was executed 2000 times with 
these dispersions. 

For this mission, because of the long (un- 
controlled) coast phase and the need for even- 
tual rendezvous by the orbiter, the apoapsis 
and inclination are the most critical param- 
eters. Figure 3 shows the final apoapsis and 
inclination for 2000 cases. For all cases the 
apoapsis is between 97 and 103 km, and the 
inclination is between 44.7 and 45.3 degrees. 

Table 2 summarizes some statistics from 
the Monte Carlo simulation. 



44 7 44 8 44.9 46 45.1 45.2 45.3 45.4 

Inclination (deg) 


Figure 3. Burnout Conditions for 2000 
cases. 


Table 2. Monte Carlo Statistics 

Variable 

Mean 

Max 

Min 


Altitude 

6.958 

7.363 

6.681 

km 

Inclination 

45.01 

45.27 

44.73 

deg 

Apoapsis 

100.22 

102.92 

97.43 

km 

Periapsis 

-3357 

-3355 

-3360 

km 

Tot. Impulse 

105.1 

106.3 

104.0 

kN*s 


5 Conclusions 

A guidance algorithm for the first stage 
of a proposed Mars Ascent Vehicle has been 
developed. This algorithm is based on a cal- 
culus of variations approach, using influence 
coefficients to drive the vehicle state to a de- 
sired terminal state. The algorithm is de- 
signed to provide good performance with very 
little on-board computation. While the exact 
configuration is subject to change, this algo- 
rithm is potentially useful across a wide range 
of applications. 

The proposed guidance algorithm has 
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been implemented and tested in a 3D0F 
Monte Carlo simulation. The results show 
that the algorithm controls the vehicle to 
relatively tight tolerances under reasonable 
environmental dispersions, keeping the final 
condition within about a quarter degree in 
inclination and three kilometers apoapsis. 
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